Density functional scheme for calculating the pair density 
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The density functional scheme for calculating the pair density is presented by means of the 
constrained-search technique. The resultant single-particle equation takes the form of the modified 
Hartree-Fock equation which contains the kinetic contribution of the exchange-correlation energy 
functional as the correlation potential. The practical form of the kinetic contribution is also proposed 
with the aid of the scaling relations of the kinetic energy functionals. 



The pair density (PD) functional theory is expected to improve upon the conventional density functional theoryQ,0, 
because the PD has a larger amount of information than the electron 

density!! H@. 

Recently the PD functional 

theory has been developed by several workers. Ziesche first proposed the PD functional theory using the natural spin 
geminals0, 0]. Another scheme was offered by Gonis with taking the two-particle densities as a basic variable^, 0|. 
Nagy generalized Gonis 's theory to the PD functional theory which has the auxiliary equation of a two-particle 
problemJlllOil. ^ rers presented the well-conceived approximate form of the kinetic energy functional within the PD 
functional theory [13. 

In this paper, we propose the other type of the PD functional theory. Our strategy is to reproduce the PD by using 
the single Slater determinant (SSD). This strategy has two kinds of merits. The first merit is that the reproduced 
PD is never unph ysic al. This is a strong merit because the 7V-reprcscntability of the PD has been an unresolved 
problem, as yetjj, LLJ, UM UM, UJ& UJ(- T ne details will be discussed later. The second merit is the interesting 
form of the single-particle equation. As shown later, it is the modified Hartree-Fock equation which additionally 
contains the kinetic contribution of the exchange-correlation energy functional as the correlation term. This form is 
reasonably intelligible to us because correlation effects are explicitly incorporated into the single-particle equation as 
the correction to the Hartree-Fock equation. Hereafter, we shall describe our scheme and the above merits in detail. 

Let us begin with the Hohenberg-Kohn theorems with taking the PD as a basic variable. The PD is defined by 

7 (2)( rr '. rr ') = ^| I JJ $+[x)$ + (rf)farf)$(x)dri&r{ |#) , (1) 

where ip(x) and ip + (x) are field operators of electrons, \I/ is the anitisymmetric wave function, and x denotes the 
coordinates including the spatial coordinate r and spin coordinate r\. We consider a system, the Hamiltonian of which 
is given by 

H = f + W + J p(r)v cxt (r)dr , (2) 

where T, W and p(r) are operators of the kinetic energy, electron-electron interaction and electron density, respectively, 
and fpyt.(r) stands for the external potential. The proof procedure is similar to the constrained-search theories |2Ctl2ll 
|ia,|23. The universal functional F[ 7 ( 2 )] is defined by 

m 7 ( 2 )i = Mm mf + wm 

= ^[7 (2) ]|r + VF|*[ 7 (2) ]) (3) 

where — ► -y( 2 ) denotes the searching over all anitisymmetric wave functions that yield a particular 7 ( - 2 * ) (rr / ; rr'). In 
the second line, the minimizing wave function is expressed as ^[7^]. Here we suppose that the minimum exists in 
Eq. ©. By virtue of the Rayleigh-Ritz principle, we getjflQ 

*U 2) ] = *o, (4) 

(2) (2) 
where ^0 is the ground-state wave function, and 7,3 (rr';rr') is the ground-state PD given by 7q (rr';rr') = 

(*o|7 (2) (rr';rr')|*o)- 
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The variational principle with respect to 7^ 2 )(rr'; rr') can be obtained by dividing the ordinary Raylcigh-Ritz 
rinciple h 

E[>yW} as 



principle into two steps: E = Min < Min C&\H 1^) >■ Using Eqs. @ and and defining the energy functional 



Eh (2) } = F[7 ( 2)] + __ / I T ( 2)(rr /. rr ') Wcxt(r)drdr ' 7 (5) 



the Rayleigh-Ritz principle is further rewritten as J 

S = Min£;M 2 )l 

7<2 > (6) 

Equations |0J and © are regarded as the extened Hohenbcrg-Kohn theorems. 

In order to perform the variational principle Eq. JBJ), it is necessary to restrict the searching area within the set 
of TV-represcntablc PD's. Only the necessary condition for the iV-rcprcscntable PD, which is called Pauli condition, 
has been realized by Coleman|ia|. However, the necessary and sufficient condition for the TV-representable PD is, 
as yet, the unresolved problem^, 0, Q, ^| 0, 0, 0, [n|. Our strategy is to restrict the searching area within 
not the set of iV-representable PD's but that of PD's which are given by SSD's. If the former and latter sets are 
respectively denoted as C and C", then CDC' holds obviously. Therefore, this strategy prevents the minimizing PD 
from being unphysical one. There is a possibility that the ground-state PD may belong to C<lC' , where C' stands for 
the complementary set of C". However, our strategy is to search the best solution within C. Best solution means the 
most appropriate reproduction of the PD with respect to searching the minimum of Eq. © within C . This spirit is 
analogous to the usual Hartree-Fock approximation because the best antisymmetric wave function is searched within 
the set of SSD's j24|. Along this spirit, we shall give the practical and useful scheme for calculating the PD. 

In order to perform the variational principle practically, we shall introduce the kinetic energy functional T s [ 7 < 2 )] 
which is defined by 

T s [ 7 < 2 >] = Min ($|f |$) , (7) 

<1> — >-y( 2 ) 

where $ is the SSD which is constructed from N orthonormal spin orbitals an( ^ ^ ~~ ¥ 7 indicates that the 

search is constrained among all SSD's which yield the prescribed 7 ( 2 '(rr'; rr'). We again suppose that the minimum 
exists in Eq. similarly to Eq. J2J). In Eq. J7J), 7 ( 2 )(rr'; rr') is expressed as the expectation value with respect to 
the SSD. Using Eq. Q), we have 

N 

7 (2) (rr /. rr / ) = _ £ / / {^( x )^( x ')^( x )^( x ') _ ipZ(x)rv(x')Mx)Mx')}dvdv'- (8) 

Taking the minimization of Eq. J7J) by the use of the Lagrange method of undetermined multipliers, the equation for 
minimizing spin orbitals can be obtained as follows: 

v \{x) + [ dx' Mr ' r0 ± Mr '' r) ^ {vft {x')^ (x'Wt (x)-p e We We (»)} 



2m 



N 



Y^s 6u tp v (x), (9) 



i/=i 



where ^(r, r') and are the Lagrange multipliers which respectively correspond to the restriction of Eq. (jSJ and 
orthonormality of spin orbitals. The Lagrange multiplier function /x(r, r') should be determined by requiring spin 
orbitals to yield a given 7^ (rr'; rr'). That is, /Lt(r, r') can be written as /x(r, r') = fj,[y^](r, r'). 
Using -F^ 2 )] and T s [7^ 2 ^], the exchange-correlation energy functional E xc [y^] is defined as 

S[ 7 (2) ] = T s [ 7 (2) ] + C/[7 (2) ] + E xc [^% (10) 

where J7[7 < - 2 ' ) ] is the usual Hartree term. Here we need the note on Eq. (|l(Jfl . The domain of F [7^] is iV-representable 
PD's , that is the set C. On the other hand, T s [7^ 2 ^] is defined with reference to PD's which belong to the set C". 
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Substituting Eq. I|10|) into Eq. JSJ) and taking the variation with respect to the PD, we can get //[7Q 2 ^](r, r') of Eq. 
©, which reproduces 7q 2 '(it'; rr'): 



SE XC [^} 



7 (2 )= ^(2) <57( 2 )(rr';rr') 



11] 



where 7q (rr';rr') is the best solution within C . 

Next we consider the exchange-correlation energy functional E xc [y^ 2 '}. Using Eq. Eq. J2J is rewritten as 



drdr 



/7 ( 2 )(rr';rr') 



(12) 



where T^ 2 )] is defined as the expectation value of T with respect to the wave function *[ 7 ( 2 )]. Note that T[^] 
coincides with the ground-state value of the real kinetic energy when 7^) = 7q 2 ' ) . Comparing Eq. 1|12[1 to Eq. 1|1U[I . 
we get 



^ C [ 7 ( 2 ^AT KC [7 (2) ]+^ C [7 (2) ], 



(13) 



where the kinetic contribution of the exchange-correlation energy functional AX^c^ 2 )] and the authentic exchange- 
correlation energy E' xc [^ 2 >\ are respectively given by 



and 



AT XC [^}^T[ 7 ^]-T S [ 7 ^} 
,2 7 ( 2 )(rr';rr') 



U[<y 



(2)1 



(14) 



(15) 



Equation l|13fl is similar to that of the conventional DFT. Here note that devising the approximate form of E' xc [y^] 
is not needed in the present scheme since E^^ 2 )] is explicitly expressed as a functional of the PD. 

Substitution of Eqs. and (|13|l into Eq. @, followed by the unitary transformation of and constant terms 
which are brought by the substitution, leads to the single-particle equation of the canonical form; 



2v72 




+ v ext (r) \xd x ) 

1 SAT^^] 



r — r' 



2 57( 2 )(rr';rr') 



1 5AT xc [j( 2 )} 



7 <2) = ^> 



2 8^f^> (r'r; r'r) 



7 (2) =7 (2) 



/V 



x {x*u{x')Xu(x')xd x )-X*u{x')xd x ')Xu{x)} = e^xd x ) 



(16) 



where Xii x ) 1S the spin orbital which is converted from ipi( x ) y i a the unitary transformation. Since the SSD is 
generally kept invariant under the unitary transformation, the reproduced form of the PD of Eq. © is rewritten as 



1 N f f 

(rr';rr') = - ^ {xfa)xt (x')x»(x)Xv( x> ) ~ X*^)xl {x')x„(x)Xv(x')} drjdr)'. 



-(2), 

7g ( rr ; rr 



(17) 



Equations (|16|l and (|17|) are final expressions for the single-particle equations which should be solved self-consistently. 
It is free from the self-interaction like the Hartree-Fock equation. 

The single-particle equation H16|) is intelligible to us because it can be regarded as the modified Hartree-Fock 
equation to which kinetic contributions are added. If AT!^^ 2 )] terms are neglected, Eq. (|16|l exactly coincides with 
the Hartree-Fock equation. Due to AT XC [7( 2 )] terms, spin orbitals deviate from solutions of the Hartree-Fock equation. 
Not only the additional terms concerning AT^^ 2 )], but also values of the Hartree and Fock potentials are different 
from the Hartree-Fock equation. All of these differences come from AT^c^ 2 '] terms, and should be recognized as 
correlation effects. 
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From another point of view, let us show again that AT XC [^ 2 ^] terms are the cause of correlation effects in the 
present scheme. Along the argument in Rcf. |22j . the following functional is introduced here: 



(18) 



where 7 < - 2 - ) (rr'; rr') is expressed as the expectation value with respect to the SSD <I>. e' [<&] can be rewritten as the 
expectation value of H with respect to $. Therefore, the minimum value of e' [$] is obtained at $ = &hf , where 
<&hf is the ground-state wave function of the Hartree-Fock approximation. On the other hand, we can also show that 
Mine' [$] = Mm {E [jW] - AT XC [7^]} similarly to the proof of Eq. (4-12) of Rcf. Therefore, we get 

{$hf\H\$ H f) = Min f^^l-AT^^ 2 ']}. (19) 
This relation indicates that the ground-state energy which is expressed as Min EVy^} in our scheme is different 

7 (2)gC" 

from that calculated in the Hartree-Fock approximation. The difference originates from AT XC [7^ 2 ']. Since AT xc [y( 2 >] 
is rewritten as Min (* T *) - Min {$ f |$) 0, it holds the inequality AT XC [7 (2) ] < 0. Therefore, it is expected 

that the ground-state energy becomes lower than ($>hf\ H \$>hf)- That is to say, correlation effects are incorporated 
into the present scheme if AT^^ 2 )] is taken into account. 

The negative value of AX^It^ 2 )] also leads to the inequality T.,^ 2 )] > T#i?[2(J, where Thf is the kinetic energy 
within the Hartree-Fock approximation. This means again that correlation effects are incorporated into the single- 
particle equation through AT XC {^ 2 '} because correlation effects have a tendency to raise the kinetic energy. In 
addition, we shall give a comment on the possibility of ATxdj^] = 0. If the wave functions which yield the PD of C 
were SSD's alone, AT XC [7^] would be always zero. However, that is not the case. Indeed, there exists the set of wave 
functions which give the PD of C and are not expressed by a SSD^|. Therefore, it is not the case that AT xc [y^] 
is always zero. 

In order to make it more feasible to perform the actual calculation of the present scheme, we propose the approximate 
form of AT xc [y( 2 >] by using the scaling relation of AT^^ 2 )] as a sum rule. The scaling relation of A T^J ^ 2 '] can be 
obtained from those of T^ 2 )] and T.,^ 2 )]. As for the scaling relation of T^ 2 '], Levy and Ziesche[27j has derived 
the relation T[y^ ] = (^ 2 T[y^]. Here 7^ (rr'; rr') denotes the expectation value of 7^' (rr'; rr') with respect to the 
scaled wave function, i.e., 7^ (rr'; rr') = C, 6 j^ (£r C r '; C r C r ')- Li the similar way to deriving the above relation, the 
scaling relation of r s [7 (2) ] is also given by T s [^ 2) ] = C 2 T„ [*y( 2 )] pif. Thus, we get AT^^ 2 '] = ( 2 AT XC [^% Along the 
technique of Ref. 29J, let li m ac t orL both sides of this relation. Utilizing the integration by parts, we obtain [3fl| 

2A?W7 (2) ]= - //^(rr'irrOr.vl^^Ljdrdr' 

7 ( 2 )(rr';rr')r.v{ / 7 ^^ (20) 

We shall use this equation as a sum rule for AT^^ 2 )]. 

With reference to the local density approximation of the conventional DFt|3|, the following form is assumed for 
AT, C [ 7 ( 2 )]: 

A7U 7 (2) ]= U At xc {^) drdr', (21) 

J J 7( 2 >=7( 2 )(rr';rr') 

where At xc (^ 2S> ) is an ordinary function of ^ 2 \ Substituting Eq. H21(l into Eq. (|2t)(l . and neglecting the surface 
integral at infinity, we get the differential equation with respect to At xc (^/^) as the necessary condition: 

4AM7^)~37 (2) 9A hy )) =0- (22) 



<9 7 ( 2) 

A similar relation for At xc (p) has been derived by Levy and Pcrdcw in the conventional DFT|3ll|. Solving the 
differential equation Q22[l. and substituting the solution into Eq. I|21|) . we finally get the approximate form of AT^^ 2 )]; 

AT, C [ 7 (2) ] = K 1 1 7 ( 2 )(rr';rr')tdrdr', (23) 
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where K is the arbitrary constant. The constant K may be available as an adjustable parameter in our scheme. 
Namely, K may be determined corresponding to the magnitude of correlation effects of the system. For example, in 
case of atomic structures, we may fix the value by fitting total energies to those of more accurate calculations j3y. 

In conclusion, this paper is aimed at searching the best solution of the PD within the set of C . The biggest merit 
of our theory is that the solution is physically reasonable. Namely, our solution is necessarily iV-representable. The 
solution even in the set C can contain the correlation effects through AT XC term more or less but definitely. This 
point is also important. It seems worthwhile to perform the actual calculation so as to check how much the present 
scheme can reflect correlation effects beyond the Hartree-Fock. This is the next step to be 

done[13- 
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